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Abstract. Inspired by a recent experiment, we study the influence of thermal fluctuations on the 
I-V characteristics of a Josephson junction, coupled to a strongly resistive environment. We obtain 
analytical results in the limit where the Josephson energy is larger than the charging energy and 
quasiparticles are absent. 



1. Introduction 

As it is well-known, the dynamics of a small Josephson tunnel junction is char- 
acterized by two energies: the Josephson coupling energy Ej and the charging 
energy Eq Q. The behavior of the junction is determined by the competition 
between these two energy scales. If the Josephson coupling energy dominates, a 
superconducting state with a well-defined phase difference across the junction is 
possible. The junction will carry a Cooper pair current I = Ij sin (j) in the absence 
of an external voltage V. Here Ij = 2eEj/h is the Josephson critical current 
(e is the electron charge). If on the other hand the charging energy dominates, 
an insulating state with a well-defined charge Q on the junction is possible. This 
gives rise to a gap in the I-V characteristics of the junction, associated to Coulomb 
blockade \ 2\: the Cooper pair current I = up to a voltage V = 2Eq / e. Here we 
defined Ec = e 2 /2C, where C is the capacitance of the junction. 

The above statements are true under appropriate conditions on the impedance 
Z{oj) of the circuit connected to the junction, as we will discuss in some detail 
below. To be specific, we will consider two configurations. The voltage-biased 
set-up, in Fig. consists of an ideal voltage source (voltage V x ), connected to 
a series arrangement of a Josephson junction and a resistance R. In a current- 
biased set-up, see Fig. an ideal current source (bias current I x ) is connected 
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to a parallel arrangement of a Josephson junction and a shunt resistance R. The 
two set-ups are equivalent if we choose I x = V x /R (3). The resistance R plays a 
crucial role and determines the detailed form of the I-V characteristics. In partic- 
ular, the way in which the transition from superconducting to insulating behavior 
of the junction manifests itself in the current-voltage characteristic depends not 
only on the parameters Ej and Ec, but also on R. 
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Figure 1. Voltage biased (a) and current-biased (b) Josephson junction. 



The voltage biased set-up, Fig.[|k, has been analyzed long ago by Ivanchenko 
and Zil'berman 01 in the limit of strong charging effects, Eq — ► oo. More re- 
cently, a perturbative approach has been developed to obtain the current-voltage 
characteristic [35]. This approach is valid for finite, but small, values of Ej / Eq ■ 
Specifically, in the zero temperature limit T = 0, a change of the position of the 
supercurrent peak is found upon increasing the ratio R/Rq, where Rq = h/Ae 2 . 
For R/Rq < 1, the supercurrent peak is centered at zero bias voltage; its width 
increases gradually as R/Rq is increased. If R/Rq S> 1, the supercurrent peak 
is centered around eV = 2Ec- This corresponds to the transition from supercon- 
ducting behavior found for small values of the external resistance to a complete 
Coulomb blockade of Cooper pair tunneling in the presence of a resistive envi- 
ronment. The range of applicability of the perturbative analysis depends on Ej, 
Ec, and R. In the low resistance limit R/Rq < 1, the approach breaks down 
at low energy (bias voltage or temperature) for any value of Ej. In the opposite 
limit R > Rq, the condition Ej/Ec <C (Rq/R) 1 ^ 2 needs to be satisfied. A non- 
perturbative analysis was developed in 0, however its validity is limited to low 
energies, eV, k B T < E c Rq/R. 

The approach discussed so far is clearly appropriate if one is interested in 
a small junction with Ej < Ec, coupled to an environment which is not too 
resistive. An example is the recent experiment by Steinbach et al. Q, in good 
agreement with the theoretical predictions of Refs. [4,6]. However, the perturba- 
tive approach is not applicable in the case of a Josephson junction with Ej > Ec, 
coupled to a strongly resistive environment, R 3> Rq, usually referred to as the 
underdamped limit Nevertheless this is a relevant situation from an exper- 
imental point of view as well. Recently, Watanabe and Haviland measured the 
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I-V characteristics of a small Josephson junction coupled to a tunable environ- 
ment [8|. In these experiments, Ej ~ Ec and the environment could be tuned 
to large values of the resistance such that R/Rq ~ 10 2 - 10 4 . No evidence of a 
supercurrent peak was observed; rather a so-called Bloch nose [9| was found for 
small values of the current and the bias voltage. This corresponds to an on-set of 
Coulomb blockade at low voltage and small current, followed by subsequent back- 
bending of the I-V curve to smaller values of the voltage at higher currents which 
indicates a cross-over to superconducting behavior. The shape of the Bloch nose 
depended on the resistance R as well as on temperature. The experimental findings 
could not be reconciled quantitatively with existing theories of Bloch nose ll9l [T0l . 
The main problem is that those theories are essentially based on the presence of 
dissipation due to quasiparticles. In the experiment [8|, the tunable environment 
is a long array of superconducting junctions in the insulating regime fTTl . and 
quasiparticles are practically absent. 

In this paper we analyze the finite temperature behavior of a Josephson junc- 
tion with Ej 3> Ec, coupled to a very resistive environment, in the absence of 
quasiparticles. We obtain, to our knowledge for the first time, explicit formulas 
for the current-voltage characteristics at finite temperature for a current-biased 
configuration with large R, see Fig. \^p. We follow the approach pioneered in 
Ref. [9], which is based on the periodicity of the Josephson coupling energy 
as a function of the phase difference cj> across the junction. As a result of this 
periodicity, the eigenstates of the junction are of the Bloch type, and the energy 
spectrum consists of bands of width of the order Ec, separated by gaps of the 
order Ej. In the limit Ej 3> Ec, the gaps are large and the junction stays in the 
lowest Bloch band at low temperature ksT <C Ej, as long as the bias current 
I x <^ Ij. If R = oo, Bloch oscillations [9] of the voltage V occur, such that its 
time-averaged value vanishes: the junction is in the superconducting state. If R 
is finite, relaxation within the lowest band modifies the dynamics of the junction. 
This results in the appearance of the Bloch nose: if I x is so small that relaxation 
prevents Bloch oscillations to occur, a finite voltage state develops with no current 
passing through the junction (Coulomb blockade), all current passes through the 
resistance. If I x is increased above a threshold value, Bloch oscillations develop 
and the a superconducting zero-voltage state is reached. Our approach is valid for 
a large temperature range and our results can be compared quantitatively with the 
experimental data of Ref. HJ. 



2. The model 



The circuit of Fig. ^> can be modelled by the Hamiltonian H = Hq + H coup \ + 
H C nv Here, 

H = ^-Ejcos4> (1) 
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describes the Josephson junction in the absence of the external circuit. The charge 
Q and the phase (/) are canonically conjugate operators, [Q, 4>] = 2ie. Since the 
Hamiltonian Q is periodic in <fi, the energy spectrum consists of bands e n (q) 
where n is the band index and q the quasimomentum associated to <p, i.e., the 
quasicharge. The spectrum e n (q) is a periodic function of q with period 2e. We 
will limit ourselves to the first Brillouin zone, — e < q < e. The corresponding 
eigenstates tp nq ((j)) are of the Bloch type. In the limit of interest here, Ej 3> Ec, 
the bands are narrow compared to the gaps. Throughout this paper we will assume 
that all relevant energies are smaller than Ej and restrict ourselves to the lowest 
energy band, n = 0, with a dispersion eo(q) = hy/EEcEj — Uq cos(-Kq/e), where 
the bandwidth is given by 



V 7T 2Eq 



(2) 



The term 



#coupl = ~TtH (3) 

couples (f) to the total current I = I x + 51, which contains the constant bias current 
I x and a fluctuating part 51. The fluctuations are induced by the environment (the 
resistance R in Fig.^i) and the dynamics of 51 is governed by H env . 

In the absence of any current, H coup i = and quasicharge q is a well defined 
variable. In the limit of small bias current I x and current fluctuations 51, the effect 
of the coupling term (|3ll on g can be analyzed perturbatively. As it was shown in 
Refs. [9 1, a quantum Langevin equation for the quasicharge operator q can be 
obtained, which reads 

dt x R dq ^ 

Here, I x plays the role of a driving force and 51 of a random force. The statistical 
properties of the latter are determined by the Hamiltonian H cnv . Specifically, for 
a resistive environment at equilibrium we have fT3l 

(6I(t)) em = , ([5I(t), 5I(t')] + ) cnViUJ = 2{hw/R) coth[hw/(2k B T)}, (5) 

where [...,.. .]+ denotes the anticommutator and (. . .) env an average with re- 
spect to H cnv . Finally, the derivative deo/dq is the voltage operator V across the 
junction, the corresponding term V/R acts as a damping term. 

We will be particularly interested in the limit of small current fluctuations. 
In this limit, quasicharge remains a well defined quantity and we can replace 
operators by classical variables in Eq. ©. In the next Section, we will analyze the 
resulting classical Langevin equation in some detail, using well-known methods 
for stochastic equations [ 12 1. We will see under which condition current fluctua- 
tions can be considered small and we will obtain the detailed, analytical form of 
I-V characteristics, as a function of the resistance R and of temperature T. 
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3. Current-voltage characteristics 

3.1. LINEARIZED L ANGEVIN EQUATION 

In the absence of fluctuations, 51 = 0, the Langevin equation (|4} has a stationary 
solution when I x = (l/R)deo/dq = V/R, with q = qo = (e/V) arcsin(I a; /i&). 
Here, = Vjj/R, with Vj, = nUo/e the maximum voltage across the junction 
corresponding to the maximum slope of eo(q). The corresponding I-V charac- 
teristic is linear, I x = V/R. In other words, all the current passes through the 
resistance R, no current flows through the junction (Coulomb blockade) and the 
phase 4> is completely undefined. This is consistent with the fact that we treat q 
as a well-defined variable. Note that the stationary solution exists only as long as 
Ix < hi f° r larger current, I x > I),, the quasicharge becomes time-dependent. We 
will discuss that case in Sec. B.2l below. 

Let us now take into account small fluctuations 51 such that the induced charge 
fluctuations 5q around the stationary solution q = qo are small. We consider the 
linear part of the average current- voltage characteristics and calculate 



Here the symbol (...) means averaging over the fluctuations of charge. In order 
to find the distribution of these fluctuations, we write q = qo + 5q, such that the 
linearized, classical Langevin equation takes the form 



where we introduced the relaxation time r such that r 1 = (ir/e)y\I 2 — J||. 
The time r defines the time scale characterizing the junction dynamics. Note that 
this is a relatively long time compared to the RC time of the circuit: for I x = 0, 
we find r ~ e/It ~ Re 2 /Uq ^> RC as Uo is exponentially small, see Eq. $2%. 

In order to solve the stochastic equation ©, we should specify the correla- 
tion function (5I(t)5I(t')}, see Eq. Q. The frequencies uj of interest here are 
small, u ~ 1/t, determined by the long time scale r characterizing the dynamics 
of quasicharge. For finite temperatures fc^T S> K/t, the correlation function 
(5I(t)8I(t'))u ~ 2k B T/R, independent of frequency and thus {5I(t)5I(t')) = 
{2ksT / R)5{t — t'). This means that the relevant current fluctuations are clas- 
sical and completely uncorrected on the long time scale r characterizing the 
quasicharge dynamics. 

The stochastic equation © can now be solved, using standard methods de- 
veloped for the analysis of Brownian motion fl2l . As a result, one obtains the 
probability distribution W(6q,t;6qo) to find the fluctuation 5q at time t given 
that 5q = 5qo at time t = 0. In the long time limit, t — > oo, this probability 



(V) = (de /dq) = V b {sm(irq/e)). 



(6) 



5q 



5 -*+5I 

T 



(7) 
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distribution does not depend on time and is independent of 5qo. It is given by the 
Gaussian distribution 



W(Sq,r) = , (8) 

V 7TJT \ 7T J 

where we introduced the parameter 7 = 2ksT/R. 

Using Eq. (J8jl we immediately obtain the average square of the charge fluc- 
tuations, ((5q) 2 } = ksTr/R. Quasicharge can be considered well-defined as 
long as these fluctuations are small, ((Sq) 2 } <C e 2 . This condition gives an addi- 
tional restriction on temperature, fe^TT /h <C R/Rq. Together with the condition 
ksT S> H/t found above we obtain the temperature window 

1 -C ksTr/h -C R/Rq. (9) 

This means in particular that R/Rq S> 1 for the analysis presented here to be 
correct. If this condition is verified, quasicharge is a well-defined variable, and 
quantum fluctuations can be neglected. For small I x the energy scale H/t is of 
the order of eVbRq/R, thus the condition © reads eVbRq/R <C ksT <C eV&. 
This has a transparent physical interpretation: the temperature has to be much 
smaller than eVf, in order not to smear the Bloch nose, but it must be larger than 
eVb(RQ I R) <C eVb in order to justify the neglect of quantum fluctuations. The 
experiment of Ref. [8 1 was performed at very large values of R/Rq, such that the 
condition © was verified. 

With the help of Eqs. © and (JSJl we readily obtain the average voltage for the 
linear part of the average I-V characteristic, 



(V) = I X R 



it k B T 



(10) 



2eR . 

Equation ( flOl represents the main result of this section. The first term on the right 
hand side of Eq. (ITOl coincides with the linear I-V curve which was obtained 
above without taking into account current fluctuations. The second term is entirely 
due to the thermal fluctuations, which suppress the resistance. In other words, the 
Coulomb blockade is smeared due to thermal activation. 

In this subsection we have considered the limit of small enough bias current, 
I x < lb, such that the Langevin equation could be linearized around a well- 
defined, stationary solution. For larger values of I x , no stationary solution exists: 
quasicharge is an oscillating function of time and the full non-linearity of the 
Langevin equation should be taken into account. In Section 13.21 we will use a 
Fokker-Planck approach to obtain the complete current-voltage characteristics. 
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3.2. FOKKER-PLANCK APPROACH 

If the applied bias current I x exceeds I b , quasicharge is a dynamical variable even 
in the absence of fluctuations. The corresponding non-linear equation of motion 
can be integrated, though, and q(t) can be found. As a result, the voltage V be- 
comes an oscillating function of time, with a period given by r. Its time-averaged 
value V over one period of the oscillation is given by [ 10 1 



V = I b R 



(11) 



In other words, on measurement time scales much longer than the period r, the 
voltage oscillations average out and a stationary situation is reached. 

In order to investigate the influence of small current fluctuations 51 on the I-V 
characteristics, the full non-linear stochastic problem © needs to be solved. This 
can be conveniently done using a Fokker-Planck approach [ 12 1. In this approach, 
the probability distribution W(q, t) to find quasicharge q at time t is found , as a 
solution of the Fokker-Planck equation 1121 



dW d 



t t ■ 

I x - I b sin I — 



dt dq VL V e / J / 2 dq 

In writing this equation, we limited ourselves again to the case of classical fluctu- 
ations only, ksT S> h/r. In addition, we note that this equation can be used only 
as long as the fluctuations are 5-correlated on the time scale defined by r. Hence, 
we need to restrict ourselves again to the temperature interval 

We are particularly interested in the stationary solution of dl 2b . W(q), reached 
in the long time limit t>r. The average voltage is then given by 

+e 

(V) = I b R J dqsm{-Kq/e)W{q), (13) 



where W{q) is normalized, W(q)dq = 1. In the long time limit t S> r, this 
corresponds to the measured, time-averaged voltage, (V) = V. 

Let us first analyze the solution of the problem in the absence of charge fluc- 
tuations [parameter 7 equals to zero in Eq. (fl2l l. Solving Eq. (fT2t and using the 
normalization condition for the distribution function we obtain 

f 5(q - q ), I x < h 

I 2e[I x /I b -sm(nq/e)] ' ±x > Jfe ' 

where q$ is the stationary solution found in Section 137X1 Substituting Eq. (fT4l into 
Eq. (fT3l for the current-voltage characteristics we immediately obtain the linear 
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Ix/Ib I/I b 




Figure 2. Calculated I-V characteristics for an underdamped junction. Left panel: bias current 
Ix/h as a function of voltage (V) /Vb across the junction, for various values of the parameter 
k B T/eV b . Curves from right to left correspond to k B T/eV b = 0, 0.004, 0.02, and 0.1. Right panel: 
current I/Ib through the junction as a function of ( V) /Vb. Curves from right to left correspond to 
the same values of ksT/eVb- 



result (V) = I X R for I < for the opposite case I x > Jj we reproduce the 
result (fTTl . The resulting I-V characteristic, a so-called Bloch nose, is shown in 
Fig.|2]as the curve corresponding to the parameter ksT/eVb = 0. 

Next we consider the current-voltage characteristics in the presence of fluctu- 
ations. To find the distribution function, W(q), we use a Fourier representation 

+oo 

W{q) = W n e i7rqn/e . (15) 

n=— oo 

Substituting the right hand side of Eq. dT3T > into Eq. ST3 and using the periodicity 
condition, W(q + 2e) = W(q), we obtain the following recurrence relation for 
the distribution function 

2 (n + i Vx ) W n = r, b (W n ^ - W n+1 ) . (16) 

Here we introduced the short hand notation tjj = 2elj/ (irj) for j = x, b. To find 
the distribution function W n in Eq. (fT6l we use the analogy with the recurrence 
relation for the modified Bessel function iTPfl 2vl v {£ > ) = — I v +\(£)]. 

Here v is the index and £ is the argument of the modified Bessel function. A 
comparison with (fT6t yields the distribution function W v -i Vx = CI v {rjb), where 
C = {2eli rjx {rn > ))^ 1 is the normalization constant. Substituting this result into 
Eq. (fT3l and using the fact that the distribution function is real, W- n = W*, we 
obtain, upon integration over q, the following result for the I-V characteristics 

e \Iir,AVb)r 

Equation dl7t is the main result of this section. It describes the complete current- 
voltage characteristics in the presence of small fluctuations. The I-V character- 
istics of Eq. ( fTTl are shown in the left panel of Fig. |2] for various values of the 
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parameter ksT/eVb- For comparison we plot the current / through the junction as 
a function of ( V) in the right panel of Fig.|2 for the same values of ksl '/eV&. This 
would correspond to a measurement in the voltage-biased set-up, see Fig.^. We 
see that the I-V curve corresponding to the limit without fluctuations is asymp- 
totically reached by the finite temperature I-V curves, both for I x / h <C 1 and 
for I x /h ^ 1- Mathematically, this is a direct result of the asymptotics of the 
modified Bessel functions [14]. Physically, this can be easily understood from 
the fact that the point I x = If, is unstable, separating a stationary quasicharge 
solution from an oscillating one. In the vicinity of lb, any small current fluctuation 
will drastically change the nature of the dynamics. Thus the I-V characteristic is 
strongly affected by fluctuations for I x ~ lb- 

We finally note that the analytical form of the result (flTt resembles the well- 
known result [15 1 for the I-V characteristic of an overdamped junction with 
R <C Rq and Ec 3> Ej. In fact, these two cases can be related by duality 
arguments. Here, we studied the temperature driven diffusion of quasicharge q, 
governed by Eq. (112b . and found the time-averaged voltage across the junction as 
the average (sin(7rg/e)) over the fluctuations of q. In an overdamped junction, it 
is the dual phase variable <p that diffuses, governed by an equation similar to (fT2b 
and one is interested in the time averaged current (sin 0) over the fluctuations of 
4>. This is why the resulting I-V curves can be related to each other, essentially 
by exchanging the role of I and V. 

4. Discussion 

In this paper we have considered the influence of thermal fluctuations on the 
current- voltage characteristics of an underdamped Josephson junction with R S> 
Rq and Ej 3> Ec, in the absence of quasiparticles. To obtain analytical results 
for the I-V characteristics, we have analyzed the appropriate Langevin equation 
using a Fokker-Planck approach. The resulting I-V characteristic is essentially a 
Bloch nose 0, smeared by thermal fluctuations. Our results are valid as long as 
the inequality © is satisfied. Since we restricted our analysis to the lowest energy 
band we need to impose in addition ksT <C Ej and I x <C Ij. 

The present work is motivated by recent experiments JH) on a small, volt- 
age biased Josephson junction coupled to a tunable environment. This environ- 
ment consisted of a long array with a large number (~ 10 2 ) Josephson junctions. 
Each junction is a small SQUID-loop, threaded by a magnetic flux Thus the 
Josephson coupling energy Ej VVSLy '($) of the array is flux-dependent and the ratio 
E^ Tay I Eq 13 * could be changed during the experiment. When <1> is tuned close to 
3>o/2, where $o is the superconducting flux quantum, the ratio E^ ia,y / E^ Tay <^ 
1 and the array is in the insulating regime. 

We ignore quasiparticles: in order for them to reach the small Josephson junc- 
tion, they have to tunnel through the entire array; the corresponding probability is 
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negligibly small. The frequency-dependent impedance Z{uj) of the array, as seen 
by the small junction, is then entirely due to the dynamics of Cooper pairs. In the 
stationary limit of interest here, we need to know the zero-frequency component 
Z(u = 0). Calculations by Efetov Q6) in the limit Ej riliy /E^ y < 1 show that 
the real part R of Z(u = 0) is strictly infinite at zero temperature; at finite T, 
thermal activation gives rise to a finite, large value of R, which depends on the 
ratio Ej Iia,y / E^ ra,y , and thus on Note that any residual quasiparticle contribu- 
tion would be insensitive to the Josephson coupling and hence independent of <£. 
Once R(&) and the parameters of the small junction are known, our results can 
be compared in principle with the experimental data of [8 1. 
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